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Abstract 

We discuss the reaction 7r“e“ —> 7r“e“7r*^ with the purpose of obtaining information on the 
77r —> vrvr anomalous amplitude We compare a full calculation at 0(jP) in chiral perturbation 

theory and various phenomenological predictions with the existing data of Amendolia et al. By 
integrating our theory results using Monte Carlo techniques we obtain a = 2.05 nb at 0{p^) and 
a = 2.17 nb after including the dominant electromagnetic correction. Both results are in good 
agreement with the experimental cross section of cr = (2.11 ± 0.47) nb. On the basis of the ChPT 
results one would extract from the the experimental cross section as amplitudes _ (g.g^i.i) 

GeV“^ and = (9.6 ± 1.1) GeV“^, respectively, which have to be compared with the low- 

energy theorem Fstt = e/(47r^F^) = 9.72 GeV“^. We emphasize the need for new data to allow 
for a comparison of experimental and theoretical distributions and to obtain .Fstt with smaller 
uncertainty. 
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I. INTRODUCTION AND OVERVIEW 


Ever since the late 1960s anomalies QiiB have played an important role in our under¬ 
standing of strong-interaction physics. Anomalies arise if the symmetries of the Lagrangian 
at the classical level are not supported by the quantized theory after renormalization, re¬ 
sulting in so-called anomalous Ward identities [3|. For the case of chiral SU(3)^ x SU(3)j:j, 
the constraints due to the anomalous Ward identities have efficiently been taken care of 
through the effective Wess-Zumino-Witten (WZW) action j^[^. The WZW action is ex¬ 
pressed in terms of the pseudoscalar octet of Goldstone bosons and contributes at 0{p‘^) in 
the momentum expansion of chiral perturbation theory (for an overview see, e.g.. Refs, 
[i, 1^). It is determined through the underlying group structure up to an overall constant 
and, in the purely strong sector, gives rise to interaction vertices involving an odd number 
of Goldstone bosons (odd-intrinsic-parity sector) j^. Using topological arguments, Witten 
has shown that the WZW action is quantized, i.e. a multiple of an integer parameter n. By 
including a coupling to electromagnetism, this parameter has been identified as the num¬ 
ber of colors Nc by comparing with the prediction of the QGD triangle diagram for the 
(anomalous) decay 7 r° —>• 77 . Once the overall factor is hxed, the (gauged) WZW action 
also predicts other interactions such as the 77 r’'' 7 r‘^ 7 r“ vertex. However, it has recently been 
pointed out by Bar and Wiese US that the W dependence in the pion-photon vertices is 
completely canceled once the W dependence of the quark charges is consistently taken into 
account. In that sense, the width of the decay 7 r° —>• 77 is predicted absolutely without ref¬ 
erence to the number of colors. The conclusion from their analysis is that one should rather 
consider three-flavor processes such as rj —>■ 7T~^7r~'j or Kj —>■ Kn to test the expected W 
dependence [lil [lll | in a low-energy reaction. However, by investigating the corresponding 
r] and rj' decays up to next-to-leading order in the framework of the combined l/N^ and 
chiral expansions, Borasoy and Lipartia have concluded that the number of colors cannot be 
determined from these decays due to the importance of sub-leadin g te rms which are needed 
to account for the experimental decay widths and photon spectra (l 2 j | . 

The decay 7 r° —77 is the prime example of an anomalous process Q and its invariant 
amplitude can be written as 


-^7r0^77 — iT ^0123 — 1- (1) 

The prediction in the chiral limit, as obtained from the WZW action, is given by 1^ Hojl 

•^.(0) = (2) 

ttFo 

where a = e^/dvr Ri 1/137, e > 0, and Fq denotes the SU(3) chiral limit of the pion-decay 
constant [^: Ft, = Fo[l + 0{mq)\ = 92.4MeV [l^. Using Eq. Q with the empirical value 
Ft, instead of Fq, one obtains for the decay rate 

r.»-.-, = ^^ = 7.73eV (3) 

in agreement with the average value of Ref. Hi: 


= (7.74 ± 0.55) eV. 


( 4 ) 


Gorrections due to explicit chiral symmetry breaking have been studied in 
E 12 , E E|- The most recent analyses yield (8.06 ± 0.02 ± 0.06) eV 17 


Refs. HE 
in SU(2)^ X 
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SU(2)^ chiral perturbation theory at 0{p^) including electromagnetic corrections at 0{e^p‘^), 
(^60±0.10)eV Q in the framework of a dispersion theory approach, and (8.10 ± 0.08) eV 
3 using U( 3 )j;^ X 11(3)3 chiral perturbation theory at 0 {p^) in combination with large- iVc 
arguments. As has been stressed in Ref. 0 , the individual experimental results show a 
large dispersion and depend on the specihc reaction used to extract the amplitude. The 
Primakoff Experiment at Jefferson Lab (PrimEx) US aims at a determination of the width 
with an accuracy of 1.5 % and will thus match the increased precision of the theoretical 
analysis. 

As mentioned above, the WZW action also predicts more complicated processes such as 
the 77 r’'' 7 r‘^ 7 r“ interaction and one clearly needs to conhrm our picture of both the leading- 
order term as well as the relevant corrections. The invariant amplitude for + 'n'~{pb) —>■ 
7 r°(p 2 ) + T^~{pz) can be written as 

(q, Pb] P 2 , Ps) = -f^37r(s2, t 2 , U 2 ; q^)e^,upae^PlP 2 Pl, (5) 


where the Mandelstam variables are dehned as S 2 = {q + Pb)"^, h = {Pb — psY, U 2 = {.'Pb — P 2 Y 
and satisfy the standard relation S 2 +U 2 +t 2 = 2M^- + MIq+q'^} The lowest-order prediction 
[0{p^)] is independent of S 2 , ^ 2 , ^ 2 , and 




9.72 GeV’. 


( 6 ) 


The physical threshold for = 0 is at = (Mj^- -l-M^o)^, +M^o), 

and <■' = M^- {Ml. - M^-M^o - Mlo)/{M^- + M^o). 

The amplitude JFstt was measured by Antipov et al. 13 Serpukhov using 40 GeV 
pions. Their study involved pion pair production by pions in the nuclear Coulomb held via 
the Primakoff reaction 

TT-+ (Z, A) ^ TT-'+ (Z, A) + 7r°, (7) 

where Z and A denote the nuclear charge and mass number, respectively. In the one-photon- 
exchange domain, Eq. m is equivalent to 


71 


7 


71 -|- 


( 8 ) 


with an almost real photon (g^ 0). Diffractive production of the two-pion hnal state is 

blocked by G-parity conservation. At CERN COMPASS j3) ^ physics program based on 


pion and kaon scattering from the nuclear Coulomb held (Primakoh scattering 3) has 
begun. The program goals include state-of-the-art measurements of the chiral anomaly 
transitions 7 r“ -|- 7 * —>■ 7r~' + and K~ -|- 7 * —>■ K~' + vr® as well as measurements of pion 
and kaon polarizabilities and radiative transitions j 3 and hybrid meson production | 3 - 
The chiral anomaly sample of Ref. [3 (roughly 200 events) covered the ranges S 2 < 10 Ml 
and 1 ^ 2 ! < 3.5 Ml < S 2 . The small ^2 range selects events predominantly associated with the 
exchange of a virtual photon, for which the target nucleus acts as a spectator. Assuming a 
constant amplitude the value 


^ Our notation and kinematics will be discussed in more detail in Sec. im 
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was extracted from the experiment 0 - The considerable discrepancy with the theoretical 
prediction of Eq. (jHI) has generated a lot of interest in clarifying the situation from both the 
experimental and theoretical sides. 

Higher-order corrections in the odd-intrinsic-parity sector of ChPT have extensively been 


studied by Bijnens et al. [15|, |27|, |28|, |29|. They included one-loop diagrams involving 
one vertex from the WZW term and tree diagrams from the anomalous 0{p^) Lagrangian 


1^. 10 , where the parameters of the Lagrangian have been estimated via vector-meson 
dominance (VMD) calculations. While the higher-order corrections to of Eq. (j21) are 
small, for they increase the lowest-order value between 7 % and 12 % within the kine¬ 
matic range of the Serphukov experiment [i^. Moreover, genuine one-loop corrections and 
0{p^) tree-level contributions were found to be comparable in size. It has also been stressed 
by Holstein [0 that the experimental value of Eq. m was obtained under the assumption of 
a constant amplitude whereas a re-analysis using a suitable dependence on the kinematical 
variables would produce a lower value [33[ 


J’g'P = (11.9 ± 0.9 ± 0.5) GeV 


-3 


( 10 ) 


and thus reduce the difference between theory and experiment. A sophisticated analysis has 
been carried out by Hannah [0 in the framework of a two-loop evaluation [0{p^)] using 
dispersive methods. From a comparison with the Antipov et al. data with the two-loop 
analysis leaving as a free parameter, Hannah obtained 

= (11.4 ± 1.3) GeVl (11) 

By also including radiative corrections, Ametller et al. |0 showed that the electromagnetic 
corrections generate a sizeable increase in the Primakov cross section, leading, in comparison 
with Eq. m, to a further decrease 

= (10.7 ± 1.2) GeVl (12) 

Finally, using an integral equation approach, Truong j0| obtained 

Jstt = 11.2GeV-^ (13) 


Further theoretical investigations of ^^ 3 ^ include calculations in the framework of dynamical 


constituent quarks |37 


The limited accuracy of the existing data in combination with the various new calculations 
clearly motivates improved and more precise experiments [ 2 ^ 1^ 10 . In a recent JLab 
experiment [ 0 , results on 7 —>■ dvr were obtained from an analysis of yp —>■ 7i~^7i^n data 
taken with the GLAS detector. The photon energy was approximately 2 GeV. A Ghew- 
Low analysis was used to extract from the cross sections over a large kinematic range. 
Preliminary results were presented by B. Asavapibhop 10 and an experimental paper is in 
preparation 
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In this present work, we will focus on the reaction 


TT + e~ 


71 + e 


■' + vr°. 


(14) 


where an incident high-energy pion scatters inelastically from a target electron in an atomic 
orbit, as shown in Fig. ^ This reaction and also K~ + e~ —>■ K~' -\- e~' + 7 r° can, in 
principle, be studied with the (190 - 300) GeV pion and kaon beams in the GERN GOMPASS 
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experiment j23|. New high-statistics pion data will allow for a determination of the form 
factor for 7 r 7 * —>■ The kaon beam three-flavor process can as well test the expected Nc 
dependence in a low-energy reaction mini- „ 

A similar (pion) experiment has already been performed at the CERN SPS ji^ . The 
experiment did not explicitly extract a value but rather claimed that the experimental 
value was consistent with theory expectations. Although the experimental backgrounds 
were described in comparisons of experimental and theoretical distributions versus 
different kinematic variables were not shown; unfortunately, the data are no longer available 
for such comparisons ji^. Without presenting such detailed comparisons, Amendolia et 
al. reported 36 events for the reaction TT~e~ —>■ 7 r“e“ 7 r° corresponding to a cross section 
of (2.11 ± 0.47) nb. However, without statistical tests such as the Kolmogorov-Smirnov 
distribution test g comparing experimental and theoretical distributions, it is not possible 
to be sure that background events were not included in the cross section value of Ref. M- 
The cross section of Ti~e~ —> 7 r“e“ 7 r°, in principle, may also include p~ production via the 
7 re —>■ pe transition However, threshold effects eliminate this background for the 300 
GeV pion beam energy of ji^, since an energy of = 600 GeV is required to produce a p 

pe. 


with rup = 770 MeV via vre 


Our work is organized as follows. In Sec. Owe briefly discuss the kinematics and formal¬ 
in Sec. uni we present the calculation of the anomaly amplitude 


ism of TT e 


TT e 7 r°. 


within SU(3) GhPT, discuss some phenomenological approaches, and use Monte Garlo 
techniques to integrate the cross section and compare it with the experimental result of 4^ . 
In Sec. HVl we summarize our results and draw some conclusions. Some technical details are 
relegated to the Appendices. 


II. KINEMATICS AND DIFFERENTIAL CROSS SECTION FOR vr e- ^ vr e-vr^ 

Following the nomenclature of Ref. |4^, the kinematics for a-|-6—>'l-|-2-|-3is shown in 
Fig. an incoming pion that scatters inelastically off an electron target and produces an 
additional 7r° in the final state: e“(pa) -|- vr“(pb) —>• e~{pi) + vr°(p 2 ) + T^~{pz)- We consider 
the target electron to be at rest and the binding energy of the electrons bound in an atom 
to be negligible relative to the incoming pion ener gy. 

We define the standard set of five invariants as 

51 = Si2 = (pi + P2f = {Pa +Pb- Psf, 

52 = S23 = iP2 + Psf = (Pa + Pb - Pl)^ 

h = tal = {Pa - Plf = {P2 + P3 - Pbf, 

t2 = tb3 = {Pb - P'if = (Pl +P2- PaY, 

S = Sab = {Pa + Pbf = (Pl +P2+ Pzf ■ (15) 

The invariant s is fixed by the incident beam energy Ei, s = -|- ml -|- 2Eime, and one is 

left with four scalar variables si, S 2 , H, and ^ 2 - 

The fourfold differential cross section, expressed in terms of the five invariants of Eq. 
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(16) 


da, is given by (see, e.g., Ref. (47j|i^ 

da 




dsids 2 dtidt 2 4(47r)‘^A(s, M 2 )(—A 4 )R 2 ’ 


where 

A(a;, ?/, z) = — 2xy — 2xz — 2yz, (17) 

and where the Gram determinant is given by 

2ml 2pa ■ Pb 2pa ■ Pi 2pa ■ Ps 
2pa-pb 2ml 2pb-pi 2pb-ps 

2pa ■ Pi 2pb ■ Pi 2mj 2pi-p3 ' 

2pa ■ P3 2pb ■ P3 2pi ■ P3 2m| 

The factor 1/16 has been extracted for later convenience. The expressions of the scalar 
products entering Eq. (HHD in terms of the invariants of Eq. m are given in Eq. (HH). 

In the one-photon exchange approximation (see Fig. El) the total invariant amplitude Af 
can be written as 

M = -iJ^3n(,S2, t2, U2] (19) 

where = eu{pi)'j^u{pa)/q'^ is the virtual photon polarization vector (g = pa — Pi) and 

Ff, = (20) 

The squared matrix element of Eq. (EnD involves the contraction of the standard lepton 
tensor known from the one-photon exchange approximation in electroproduction processes, 

W = {^PaPl + 2piPa + q^g^'") , (21) 

with the hadronic tensor and is given by 

\^3n\‘'W^F,F,= \F3^\‘^[Apa-Fpi-F + q^F-F). ( 22 ) 

The explicit expression for Eq. (j22|l is given in App. ^ 



III. THEORETICAL DESCRIPTION OF THE AMPLITUDE 

In this section we describe the theoretical input to our analysis of the reaction n~e~ —> 
7r“e“7r°. Since we work in the one-photon exchange approximation, it is sufficient to consider 
the transition-current matrix element 

(7r°(p2),7r"(p3)|Jp(0)|7r“(pfe)) = A3^(S2, ^2, M2; g^)e/.i/p,7PbP2P3i 


^ Our normalization of the electron states and of the Dirac spinors is given by 


{p', s\F s) = 2E{p){2Tr)^6^{p' - p)Ss's, u{p, s)u{p, s) = 2me8s's- 
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where is the electromagnetic current operator (including the elementary charge). In the 
isospin-symmetric limit, ^^ 3 ^ is a completely symmetric function of the Mandelstam variables 
S 2 , t 2 , and U 2 . In the physical region, the Mandelstam variables satisfy the standard relation 
S 2 + U 2 + t 2 = 2M^_ -|- M^o + We will lay emphasis on a calculation within the framework 
of chiral perturbation theory at 0{p^) but will also discuss the results of some (more) 
phenomenological approaches. This will allow us to have an estimate of effects which would 
be subsumed in terms of 0{p^) and higher. 


A. Chiral perturbation theory at 0{p^) 

Besides the neutral-pion decay into two photons, the amplitude for y-t-vr” —>• 7 r°-|- 7 r~ is of 
prime interest for testing our understanding of anomalous Ward identities. In the momentum 
and quark-mass expansion, its leading-order contribution is of 0{p'^) and originates from the 
Wess-Zumino-Witten action j^|^. The interaction Lagrangian relevant in the presence of 
external electromagnetic helds (described by the vector potential M^) is given by [ 1 , |^, 

•Cwzw = -eAp" + !^€'“''”S„.4,.4,TVl2l32(ra„f/t - Wd^U) - QWQd^U + QUQa^U% 

(23) 

where Q = diag(2/3, —1/3, —1/3) denotes the quark-charge matrix and 

U = exp ) </ = I] ^a(pa, (24) 

contains the Goldstone boson helds. The current 

=- -Tr(QdMU^dpUU^d^UU^ + QU^dMU^dpUU^dM), (25) 

loTr"* 

by itself is not gauge invariant and the additional terms of Eq. (P?|) are required to obtain 
a gauge-invariant action. The hrst term of Eq. (0!^ gives rise to the 3(/) -|- 7 coupling (see 

Fig. El) 

^wzw = g^ 2^3 {Qd,,(j)dp(j)d^(j )), (26) 

whereas the second is responsible for the 7 r° —77 decay not discussed in this paper. 

Weinberg’s power counting scheme establishes a connection between the chiral ex¬ 
pansion and the loop expansion. Since the anomalous sector only starts at 0{p‘^), the con¬ 
tribution at 0{p^) results from either one-loop diagrams with exactly one vertex from the 
WZW term or tree-level diagrams with exactly one vertex from the anomalous Lagrangian 
at 0{p^). In order to determine the one-loop contributions we need, besides Eq. the 
WZW contribution involving 5 Goldstone bosons, 

^wzw = (27) 

and the lowest-order Lagrangian in the presence of an electromagnetic held, 

£2 = + up), ( 28 ) 
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where the relevant covariant derivative is given by D^U = d^U + U] and y = 2BqM 

contains the qnark-mass matrix M and Bo is related to the scalar qnark condensate in the 
chiral limit. The most general anomalons Lagrangian at 0{p^) has recently been derived 
in Refs. According to Table V of Ref. [3l|, seven strnctnres have the potential of 

contribnting to 30+7 vertices. In principle, the corresponding low-energy conpling constants 
shonld be calcnlable from the nnderlying theory. However, since we cannot yet solve QCD, 
the parameters are either taken as free parameters that are htted to experimental data or 
are estimated from models snch as meson-resonance satnration |4^ l5nl|. 

In what follows, we will make nse of the SU(3) version of chiral pertnrbation theory, 
becanse this will allow ns in fntnre calcnlations to make contact with other anomalons 
processes involving in addition kaons. Moreover, we note that previons calcnlations at 
0{p^) were performed for real photons, = 0, becanse the amplitnde was embedded 
in a Primakoff reaction, where the virtnal photon of the Conlomb held of a heavy nnclens 
is qnasi real. For onr reaction snch an approximation is not admissable which we will also 
explicitly see when we discnss the results. 

The relevant one-loop diagrams are shown in Figs. El El and El and fall into two distinct 
groups. The hrst category just includes one graph (Fig. E)) whose loop is attached to one 
single vertex (the loop is composed by one internal line so to speak), while in the second 
category the loop always binds two different vertices together (the loop is therefore composed 
of two internal lines). Moreover, at 0{p^) one obtains a contact contribution shown in Fig. [71 
Combining the 0{p‘^) and 0{p^) results, multiplying with a factor of for each external 
pion line, and renormalizing the coefficients of the O(p^) Lagrangian, the result of the one- 
loop calculation in SU(3) ChPT at 0{p^) is given by 


^3n{s2,t2,U2]q‘^) — 


47r2F3 


1 + Cm^MI + Cq2q^ + 


1 I ^2 + ^2 + ^2 ( h 

327r2F2 1 3 I m 


+q^ In 
4 r 


+ g(s2 + U2 + t2 + 3q^) 


3 f 


F(s2 , M^^) + F{t2, M^) + F{u 2, Mi) + 3F(g0 M^) 


(29) 


The constants Cmi and Cq 2 are linear combinations of renormalized low-energy coupling 
constants of the most general odd-intrinsic-parity Lagrangian at 0{p^) HHIHIil, 


Cm2 


Cq2 


512Tr‘^{Lli - L 


6,e 

14 



51271^ 

3 


(C 


6,e 

13 




(30) 


These coefficients still depend on the renormalization scale p but in such a way that the 
scale dependence of the logarithms in Eq. (Ei is precisely compensated. The function F 
originates from a standard one-loop integral of mesonic chiral perturbation theory and is 
given by 

F{a, m^) 

J^°\x) 
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= 


t(0) (_ « 

\m^) 2 ’ 

= [ dzln[l + — z) — iO^] 

Jo 










r 2-aln(^) (x<0) 

= < -2 + 2 ^I - 1 arccot - l) (0 < x < 4) , 
[ -2-aln(^)-ivra (x > 4) 


with cr(x) = — 4/x for x ^ [0,4]. 

Up until now, we have carried out everything which is necessary to meet the requirements 
of a consistent 0{p^) calculation within the framework of mesonic ChPT. However, our result 
still depends on two unknown parameters {Cmi and Cq 2 ) which prevent us from predicting 
observables such as the total cross section or distributions. Of course, these low-energy 
coupling constants (LECs) can in principle be determined within appropriate experiments. 
Here, we actually move on and estimate the constants by using theoretical means which 
evidently have to go beyond mesonic ChPT. The LECs are supposed to include whatever 
QCD information on all particles which do not belong to the Goldstone boson octet. At low 
energies the lightest are expected to be significant and we are thus led to consider the effects 
due to the vector mesons |5^. For that purpose we made use of the nonlinear chiral 
Lagrangian of Ref. js^, evaluated the tree-level diagrams contributing to 7 * -|- 7 r“ —-|- 7 r“ 
involving internal vector meson lines, and expanded the propagators to be able to collect 
the arising 0{p^) pieces (for later purposes we also consider the expressions keeping the full 
propagators). A comparison (matching) with the polynomial p® pieces of the most general 
anomalous Lagrangian yields the estimate 



3 

1024:7i‘^mY 


T 6,e T 6,e t 6,e q r 6,e 

-^13 ? -^6 — “-^14 — ^-^5 • 


Using mv = rnp = niuj in SU(3) results in 


Cm2 



(31) 


(32) 


Instead of expanding the vector-meson propagators in the vector-meson saturation calcula¬ 
tion, we could also keep the complete propagators. This would correspond to the replacement 


Cm^M^ + Cq2q^ 


3M^ + 4g^ -^2 G U2 \ 3 

2m2 2 - S 2 m'j, - t2 - ^ 3 / 2 


(33) 


in Eq. We then obtain some estimate of higher-order terms beyond 0{p^). 

A full calculation of all C(e^) radiative corrections as well as the isospin symmetry break¬ 
ing effects due to the different u- and d-quark masses is beyond the scope of the present 
paper. As discussed in Ref. jssf, the most important electromagnetic correction originates 
from a photon-photon fusion into a neutral pion (see Fig. E> yielding an additional contri¬ 
bution of the type 




2e^Fi 


4n-2FJ 


(34) 


Due to the 1/^2 pole of the exchanged photon, Eq. (j34|l becomes important for small values 
of ^ 2 - Note that including the single diagram of Fig. |S1 does not lead to a conflict with gauge 
invariance. On the other hand, other electromagnetic corrections were found to be very 
small in Ref. [s^, where the authors concluded that their full calculation can be very well 
reproduced by adding only the contribution of Eq. (El- 
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B. Phenomenological approaches 


In our analysis we will also compare data with phenomenological calculations using three 
different forms of the transition-current matrix element: 

1. Phenomenological ansatz of Terent’ev j^: 


J^3n{s2,t2,U2;q‘^) = 


1 -|- A, 


S2 


+ 


- S2 




h 


+ 


U2 


mt 


U2, 


+ Ai, 




(35) 

Here, the multiplicative constant refers to the low-energy prediction of Eq. © 
and the variation of the function is supposed to come mainly from vector-meson- 
exchange diagrams. The parameters Ap (A^) implicitly contain factors of and 

are related to the partial widths r(p+ —> and r(p — >• vry) [r(a; —>• e+e“) and 

r(a; —>■ Svr)]. In the numerical analysis we make use of rup = 770 MeV and = 782 

MeV. Comparing Eq. (pH|l with Eq. (p?|l would result in Ap 1/2 and A^ 3/2 in 
comparison with Ap < 1/2 and A^ 2.6 of Ref. j2lj. 


2. Pole model of Rudaz including vector-meson dominance (VMD) j42| |: 


( 0 ) 


J^3As2,t2,U2;q ) = 


m,. 


m~ 


- S2 


+ 


mj^ -12 


+ 


m^p - U2 


mlj — q^ 2) 

3. VMD model including the effects of final state p-wave vrTr scattering 0: 


(36) 


_ , . 1 e 

J'3Tr{S2,t2,U2) = —- 


X 


2 47r2F3 

TT 

1 — S2/m 


1 - 


Di{s2) 


mz 


m^p - S2 


+ 


mz 


ml - t2 


+ 


mz 


ml-U2^ 


l-t2/ml'^ fl-U2/m^ 

Di{t2) I i D,(u2)‘ 


(37) 


where 


Di(a) = 1-^ 


m2 967r2i72 


In 


m; 

M2 


247r2F2 


F(a,MZ). 


(38) 


Note that the ansatz of Eq. (jSZD has only been derived for real photons, q^ = 0, and we 
therefore have to expect some shortcomings in the description of 7 * -f 7r“ —>■ tt^ + 7r~. 


C. Results and discussion 

Using Monte Carlo techniques we determined the total cross section based on Eq. 
for kinematical variables inside the region specified by 


0.0184 GeV^ < Si < 0.186 GeV^ (39) 

0.0754 GeV^ < S 2 < 0.325 GeV^ (40) 

-0.236 GeV^ < ti < -0.001 GeV^ = (41) 

-0.269 GeV^ < G < 0, (42) 
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which are the minimal (maximal) values obtained from the equations for the kinematical 
boundaries [see Ref. 0 and Eq. dCl]. For the generated invariants the positivity of —A 4 
of Eq. (fTT)]) is checked and events with positive A 4 are rejected. In order to check the Monte 
Carlo calculations we also compared the result with an explicit numerical integration using 
the simplihcations of a constant T^-k and 0 (see Appendix [01) . 

The results for the total cross section are shown in Table Q] The hrst column denotes 
the model/theory and the corresponding parameters used; the second column contains the 
integrated cross section for each case with hxed to e/(dvr^Fj) = 9.72 GeV“^. The cases 
1), 4), 5), and 6 ) of Table HI were already used in Fig. 4 of Ref. and our corresponding 
cross sections are in reasonably good agreement. In the third column we show the respective 
physical threshold amplitudes. In the chiral limit, the threshold amplitudes should reduce to 
the low-energy prediction e/(dvr^Fg ) of Eq. (jH)). In this context it is important to realize that, 
in general, the dependence of the threshold amplitude on results from both kinematical 
variables and explicit symmetry breaking. The models of 2) - 6 ) can only account for the 
hrst type of dependence and the corresponding modihcation is of the type const, x 
where the relevant constant depends on the parameters of the model and is of the order of 
1 . The results of 7) - 11), in addition, contain corrections from Goldstone boson loops of the 
form const.' x /(dvrE^r)^ which are of both kinematical and chiral symmetry breaking type. 

In the calculation of 9), the l/G singularity due to the electromagnetic correction of Eq. 
(El generates a 17 % increase in units of e/( 47 r^F,^), but when integrated over G ultimately 
leads to a less pronounced contribution to the cross section. Finally, in the last column we 
have also included the overall factor which one extracts based on the experimental 


result (2.11 ± 0.47) nb of Ref. 4^ if one treats Fg))'’ as a free parameter in the respective 


model. The error in only reflects the error in the experimental cross section and 

does not include any error estimate implied by the models. However, we explicitly do not 
advocate such an extraction as a strict test of the low-energy theorem of Eq. (0, because 
it introduces a bias in how the chiral limit is approached. Rather, at this point, the only 
rigorous approach consists of using the chiral expansion as in Eq. Ei and confronting it 
directly with experimental results. 

Using the estimate of Eq. El for the parameters Cm^ and Cq 2 , we obtain as the GhPT 
result at 0{p^) 


a = 2.05 nb. 

By also including the most prominent electromagnetic correction |3 
photon fusion into a neutral pion the result increases slightly: 


(43) 

in terms of photon- 


(T = 2.17 nb. 


(44) 


Both results are in excellent agreement with the experimental result (2.11 ± 0.47) nb of 
Ref. 113 . In general, the conceptual advantage of the GhPT calculation over the remain¬ 


ing empirical models is that it is the only calculation which naturally incorporates both 
genuine quantum effects (loops) and higher-order corrections (as estimated from the VMD 
saturation). Moreover, in principle, a controlled improvement is possible by performing a 
complete 0{p^) calculation, whereas the remaining calculations suffer from the absence of a 
systematic method of improvement. A comparison of 8) with 7) clearly shows the necessity 
to include the consequences resulting from the virtuality of the exchanged photon. This can 
also be seen from the transitions 2) to 3) or 4) to 5) in the calculations using Terent’ev’s 
model. On the other hand, the spread of the obtained cross sections is an indication that 
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higher-order terms (in the chiral expansion) may play an important role. This conjecture 
is supported by an analysis for a pion beam energy of 150 GeV which leads to total cross 
section results between 0.17 nb and 0.21 nb, i.e., the results scatter substantially less than 
for the higher energy. 

In Fig. ini we show the generated distributions of events as functions of the invariants Si, 
S 2 , ti, and t 2 as obtained using chiral perturbation theory at 0{p^) [see Eq. (j^HI)] with the 
low-energy constants of Eq. (El and the most prominent electromagnetic correction of Eq. 
El- The results are based on the generation of 100 000 events restricted to the kinematic 
region 


2Ml <si< lOM^, 

AMl <S2< lOM^, 

-0.015 GeV^ < G < -0.001 GeV^ 

-0.269 GeV^ < ta < 0. (45) 

The regions have been chosen such as to avoid the 1/ti pole and to be far enough away from 
the p-meson production threshold. In Fig. we show the differential cross section da/dti 
as a function of ti for three different choices of the low-energy constants of Eq. (El. Glearly, 
My = 0.2 GeV (triangles) would correspond to unrealistically large higher-order terms which 
is supported by the drastically different behavior for this case. Figure ITTI illustrates how the 
inclusion of the electromagnetic correction of Eq. El affects the differential cross section 
da/dti- The calculational errors in the plotted points of Figs. HU] and [TT] range from 0.1 % at 
ti = —0.015 GeV^ to 0.3 % at = —0.001 GeV^. Finally, Fig. [T21 contains a comparison of 
different (model) calculations for the differential cross section da/dti. Here, the calculational 
errors in the plotted points range from 0.03 % at H = —0.015 GeV^ to 0.07 % at = —0.001 
GeV 2 . 


IV. CONCLUSION 


We have studied the reaction 7r~e~ —>■ 7r~e~7r^ with the purpose of obtaining information 
on the yvr —tttt anomalous amplitude ^^ 3 ^. In Tabled we have summarized the results of 
various phenomenological models and of a full calculation at 0{p^) in chiral perturbation 
theory for the total cross section using the kinematical conditions of Ref. 4^1 . In particular, 
by integrating the GhPT results using Monte Carlo techniques we obtain a = 2.05 nb at 
0{p^) and a = 2.17 nb after including the dominant electromagnetic correction. Both 
results are in good agreement with the experimental cross section of cr = (2.11 ± 0.47) 
nb li^. On the basis of the ChPT results one would extract from the the experimental 
cross section as amplitudes = (9.9 ± 1.1) GeV“^ and ^^ 3 °^'^^*'^ = (9.6 ± 1.1) GeV“^, 

respectively, which have to be compared with the low-energy theorem = e/( 47 r^F(f) = 
9.72 GeV~^. We emphasize the need for new data to allow comparison of experimental and 
theoretical distributions and to obtain with smaller uncertainty. In order to further 
support our Endings and to obtain with a smaller uncertainty, it would be useful for 
future experiments to also consider distributions such as those shown in Figs. IHl-IT^ 
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APPENDIX A: SCALAR PRODUCTS 


The ten scalar products appearing in the calculation of the differential cross section of 
Eq. (Ilf)|l may be expressed in terms of the five standard invariants of Eq. (HH) as 0 


2 Pa ■ Pfe = s-ml-ml, 

2 Pa-Pi = ml + ml-ti, 
2pa-P2 = si + ti 
2 pa-P 3 = S-Si+t 2 -ml, 
2 pb-pi = s-S 2 + ti- ml, 


‘ 2 pb-P 2 = S2 + t2-ti- ml, 

2 pb-p 3 = ml + ml-t 2 , 

2pi-p2 = Si—ml — ml, (Al) 

2 Pi ■ p3 = s- Si-S 2 + ml, 

2 P2 ■ P3 = S2-ml- ml- 


APPENDIX B: THE SQUARED AMPLITUDE 


In order to evaluate Eq. ( 1221 ), one makes use of the double epsilon expression 


p'cr' 


—det(p"“ ), a = fi,u,p,a, a'= p',v', p',a'. 


The evaluation of the expression in brackets in Eq. ( 1221 ) is straightforward but tedious and 
we only quote the hnal result: 

(^Pa- Fpi ■ F + ]^q^F ■ F^ 

= (^Pa ■ Pi + [Pb(p2 ■ P3f + P2(P6 ■ P3f + pl{Pb ' P2)^ “ pIpIpI “ ‘2pb ' P2Pb ' P3P2 ' Ps] 

+ 2 pa ■ PbPl ■ Pb[plpl - (P 2 ■ P 3 )^] 

+2p„ ■ P 2 P 1 ■ P2\pIpI - {pb ■ P3f] 

+2pa ■ PsPl ■ P3\pIpI - {Pb ■ P2)^] 

+2(p„ • PbPi ■P2 +Pa- P 2 P 1 ■ Pb) {Pb ■ P 3 P 2 ■P3-Pb- P2pI) 

+ 2 {pa ■ PbPl ■P 3 +Pa- P3P1 ■ Pb) {Pb ■ P2P2 -Ps-Pb- P3P2) 

+ 2 (Pa ■ P2P1 ■ P 3 + Pa ■ P3P1 ■ P 2 )(P 6 ' P 2 Pb ' PS “ P 2 ' P 3 Pb) ■ (Bl) 


We deliberately did not express the scalar products pf, pi, and p| in terms of the masses, 
because it is then straightforward to apply Eq. jmi for other processes involving particles 
with different masses, such as, e.g., 7 *(g) + K~{pb) —> 7r°(p2) + ^'“(ps). Note that Eq. (jBljl 
is manifestly symmetric under both the exchange pi •«-*• Pa and the exchange of any two 
elements of {P 6 ,P 2 ,P 3 }- Finally, using Eqs. (lAlf) the scalar products in Eq. ED may be 
expressed in terms of the invariant variables of Eq. m- 
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APPENDIX C: DIRECT CALCULATION 


Using the covariant normalization for both bosons and fermions, the differential cross 
section for e~{pa) + '^~{Pb) —e“(pi) + vr°(p 2 ) + can be written as 


da = 


\liPa-PhY -mlMl (27r)® 


^ -P2(Cl) 


2Ei 2E^ 2Ei ’ 


where we consider the final-state pions in their rest frame (denoted by *) and the ejected 
electron in the laboratory frame. Integration with respect to pg and E 2 yields 


da = 


4 Y^(pa -PbY - mlMl (27r)'’ 
Let the pion beam define the positive z axis. Using 

d?pi 


m 


■d^Pl |P2M^2 

2Ei 


/o El 


= irdEidpiz 


in combination with 


2|pfe|pi^ = S 2 -MI- + s- Ml), 


2meEi = 2m ^ — ti, 

me\Pp\ = \J{Pa-Pb)'^ -mlMl, 


we obtain 

and can thus write 
da = 


ndEidpiz = 


TT 


dtids2 


4 iJ{pa-Pb)'^ -mlMl 

1 I jii^ I 7rdt\ds2 


\P2\d^l 


^\J{Pa- Pb)"^ - rnlMl (27r)'’ ^\J{pa ■ Pb)"^ - mlMl 


(C2) 


To obtain the kinematic boundaries of S 2 and ti one can either solve for ti in terms of S 2 , 
or vice versa ji^ . 

= ml + mj- ^[(s + ml- ml){s - S 2 + mj) T AC2 (s, ml, ml)X^/^{s, S 2 , mj)]. 


s- = s+ m~, - ^^[{s+ ml-ml){ml + ml-ti)T >^^^‘^{s,ml,ml)X^^‘^{ti,ml,ml)], 


’2 


2ml 


(C3) 


where X{x,y,z) is defined in Eq. ffTTjl . 

In order to test our numerical integration programs, we evaluated Eq. (Q under the 
following simplifying assumptions: We neglected terms containing the square of the electron 
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mass and we assumed T'i-u to be constant [see Eq. ©]• Under these assumptions we obtain 
for the angular integral 


J d^KApa-Fpi- F + tiF ■ F) = ~ti{s2-AM^)[tj + 2ti{s-S2- 
Inserting Eq. (IDl into Eq. (ini and using Eq. (I22D results in 


M^^)+{s-M^^r+{s-s2n 

(C4) 


da 


dtids 2 (s — M‘^)‘ 


■\eF, 


Stt I 




y/^ti 


[t^+ 2 ti{s—S 2 —M^) + {s—M^) +(s—S 2 ) ] 


(C5) 

for ml ^ 0 and Fsn = const. We hrst integrate Eq. (insi) with respect to ti from = 
— (s — M^){s — S 2 )/s to ti = where is an experimental cut on the maximal value of tp. 


rp, = 1 1 |,^ |2 fa-4M;p 

Jt- ^ dtids2 [s — 

X - tf) + 2(s - S 2 - - tl) + [(s - MlY + (s - S 2 )^] In 

(C6) 

where AM^ < S 2 < S 2 = s[l + 1 \/(s — M^)] as implied by the physical boundary conditions. 

A numerical integration of 



(C7) 


yields a{tl) = 1.864 nb for = —0.001 GeV^ which agrees within less than 0.5 % with the 
Monte Carlo result 1.855 nb (we used = 139.57 MeV). 
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FIG. 2: One-photon exchange approximation of the reaction e {pa) +t^ {pb) —> e (pi) + 'k^{p 2 ) + 
7I'“(P3)- 



FIG. 3: WZW diagram obtained from Eq. (|^ . 



FIG. 4: One-loop diagram obtained from expanding the first term of Eq. (ESI) to fifth order in the 
Goldstone boson fields and contracting two lines to form a loop. 
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FIG. 5: One-loop diagram obtained from contracting two lines of Eq. (EH) with the two lines of 
the 2 (/> -|- 7 vertex from £2 to form a loop. 



FIG. 6 : One-loop diagrams obtained from contracting two lines of Eq. (EH) with two lines of the 
A4> vertex from £2 to form a loop (cuts in the S 2 , t 2 , and U 2 channels, respectively). 



FIG. 7: Contact diagram obtained from £ 5 . 
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FIG. 8: Electromagnetic correction. 


TABLE I: Total cross section for 7r“e“ —> 7r“e“7r*^ as obtained in different models and chiral 
perturbation theory (see text). The third column contains the physical threshold amplitudes for 
= 0. In the last column we made use of the experimental result (2.11 ± 0.47) nb of Ref. 
to convert this into the value for which would be extracted based on the given model: 

= 9.72 ^ cJexp/o'modei GeV“^. Here we follow the common practice to quote the extracted 
values. The error only reflects the experimental error and does not include any error estimate 
implied by the models. We stress that in the framework of ChPT the overall factor of Eq. (EH) is 
not a free parameter. Therefore, the “extracted” values for the ChPT calculation have to be taken 
with a grain of salt (see the discussion in the text). 


Model/theory 

cross section 

.^(O)extr 

Stt 


[nb] 

[GeV-3] 

[GeV-3] 

l)-^3. = 4;;% = 9.72GeV-3 

1.92 

9.7 

10.2 ± 1.1 

2) Terent’ev, Eq. with Ap = 0.5 and = 0 

2.80 

10.3 

8.4 ±0.9 

3) Terent’ev, Eq. with Ap = 0.5 and A^^ = 1.5 

2.62 

10.3 

8.7 ± 1.0 

4) Terent’ev, Eq. with Ap = 0.35 and A^ = 0 

2.51 

10.1 

8.9 ± 1.0 

5) Terent’ev, Eq. ((HH with Ap = 0.35 and A^ = 3.2 

2.18 

10.1 

9.6 ± 1.1 

6 ) Rudaz, Eq. (IHH|) 

2.36 

10.0 

9.2 ± 1.0 

7) ChPT at C>(p®) [Eq. (EH] without q‘^ dependence 

2.33 

10.4 

9.2 ± 1.0 

8 ) ChPT at 0{p^) [Eq. (EH] ''^ith q‘^ dependence 

2.05 

10.4 

9.9 ± 1.1 

9) ChPT at 0{p^) [Eq. (EH] ''^ith q‘^ dependence 




plus electromagnetic correction of Eq. (jHll) 

2.17 

12.1 

9.6 ± 1.1 

10) ChPT at 0{p^) with modified dependence of Eq. (jH^l) 

2.83 

10.5 

8.4 ±0.9 

11) Holstein, Eq. (|37|) 

3.05 

10.4 

8.1 ±0.9 
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FIG. 9: The distributions of events as functions of the invariants (a) si, (b) S 2 , (c) ti, and (d) t 2 
as obtained using chiral perturbation theory at 0{p^) [see Eq. (j29l) ] with the low-energy constants 
of Eq. (I, ‘121) and including the most prominent electromagnetic correction of Eq. (El). 
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FIG. 10: Differential cross section do/dt\ as a function of ti using chiral perturbation theory at 
0{p^) [see Eq. (|^ ] including the most prominent electromagnetic correction of Eq. (jni. The 
low-energy constants of Eq. (El have been fixed using My = 0.2 GeV (triangles), My = rup 
(circles), and My = 2mp (squares), respectively. 



EIG. 11: Differential cross section do/dti as a function of ti using chiral perturbation theory 
at 0{p^) [see Eq. ([21] with the low-energy constants of Eq. El without (circles) and including 
(triangles) the most prominent electromagnetic correction of Eq. (1(141) . 
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FIG. 12; Differential cross section da/dti as a function of ti for the model of Terent’ev, Eq. (l.‘-{5jl 
with Ap = 0.35 and = 3.2 (triangles), ChPT at 0{p^) [Eq. ((201)] plus electromagnetic correction 
of Eq. (ini|) (circles), and the model of Holstein, Eq. (1^71) (squares). 
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